Symmetry-induced interference effects in metalloporphyrin wires 
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Organo-metallic molecular structures where a single metallic atom is embedded in the organic backbone are 
ideal systems to study the effect of strong correlations on their electronic structure. In this work we calculate 
the electronic and transport properties of a series of metalloporphyrin molecules sandwiched by gold electrodes 
using a combination of density functional theory and scattering theory. The impact of strong correlations at 
the central metallic atom is gauged by comparing our results obtained using conventional DFT and DFT+U 
approaches. The zero bias transport properties may or may not show spin-filtering behavior, depending on the 
nature of the d state closest to the Fermi energy. The type of d state depends on the metallic atom and gives rise 
to interference effects that produce different Fano features. The inclusion of the U term opens a gap between 
the d states and changes qualitatively the conductance and spin-filtering behavior in some of the molecules. 
We explain the origin of the quantum interference effects found as due to the symmetry-dependent coupling 
between the d states and other molecular orbitals and propose the use of these systems as nanoscale chemical 
sensors. We also demonstrate that an adequate treatment of strong correlations is really necessary to correctly 
describe the transport properties of metalloporphyrins and similar molecular magnets. 

PACS numbers: 31.15.A-, 73.63.-b, 72.80.Ga 



I. INTRODUCTION 

A key issue in the field of molecular electronics is the 
search for molecular compounds that give rise to new or im- 
proved functionalities. Porphyrin molecules constitute in this 
respect promising candidates and are as such receiving in- 
creased attention. These molecules play an essential role 
in many biological processes such as electron transfer, oxy- 
gen transport, photosynthetic processes and catalytic substrate 
oxidation. 1 . Porphyrins have been extensively studied in the 
past by biologists and chemists 2 ^—. However, an increas- 
ing number of theoretical 6-8 and experimental 910 physics an- 
alyzes have appeared in the past few years. Progress in 
the design of supramolecular structures involving porphyrin 
molecules has been rather spectacular—. 

Porphyrins are cyclic conjugate molecules. Their par- 
ent form is the porphin, which is made of four pyrrole 
groups joined by carbon bridges and has a nearly planar D4h 
■ symmetry^ 2 -. The alternating single and double bonds in its 
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structure, make this molecule chemically very stable. In por- 
phyrin systems, the porphine is the base molecule, but differ- 
ent functional groups can join the macrocycle by replacing a 
peripheral hydrogen. In addition, the macrocycle can accom- 
modate inside a metallic atom, such as Fe (which is the base 
of the hemoglobin in mammals), Cu (hemolymph in inverte- 
brates), Mg (chlorophyll in plants), etc. Large porphyrin sys- 
tems can undergo certain ruffling distortions because of its pe- 
ripheral ligands, its metallic atom inside or the environment 13 . 
It is precisely the large number of possible configurations, 
which give rise to a wide variety of interesting properties that 
make these molecules very attractive for molecular-scale tech- 
nological applications. 

In this work, we present a theoretical study of the elec- 
tronic and transport properties of porphyrins sandwiched by 
gold electrodes using Density Functional Theory (DFT). The 



porphyrine molecule is attached to the gold surface by a thiol 
group and oriented perpendicular to it, as sketched in Fig. {T). 
The metallic atom placed inside can be Fe, Co, Ni, Cu and Zn. 
For the sake of comparison, we have also studied the porphine 
compound, which has no metallic atom. Since DFT fails to 
describe correctly transition metal elements in correlated con- 
figurations, we have also adopted a DFT+f/ approach. We 
show that the inclusion of strong correlations is necessary to 
accurately characterize the transport properties of some of the 
metalloporphyrin complexes. 

The outline of the article is as follows: in section II we 
briefly explain the DFT+U flavor used in our calculations. In 
section III we give the computational details and in the last 
two sections, IV and V, we present and discuss the electronic 
and transport properties of the junctions, respectively. 



II. DETAILS OF THE DFT+ U APPROACH 

DFT has emerged as the tool of choice for the simula- 
tion of a wide array of materials and nanostructures. How- 
ever, the theory fails to describe strongly correlated electron 
systems, such as embedded 3d transition metal or 4f rare 
earth elements. Apart from the fact that even exact DFT 
can not describe all excited states 14 , the approximations in- 
cluded in the exchange-correlation potential induce qualitative 
errors in correlated systems. There have been many attempts 
to fix these problems. These include generating exchange- 
correlation functionals specifically tailored to the system un- 
der investigation, but then those are frequently not transfer- 
able. Another attempt is based on removing the electronic 
self interactions introduced by the approximate treatment of 
exchange in DFT 15 . This unphysical self-interaction is a sig- 
nificant source of errors when electrons are localized. How- 
ever it is not clear that just by removing the self-interaction 
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FIG. 1 : (Color online) Schematic view of a Fe-porphyrine molecule 
between gold leads. Yellow, dark yellow, grey, black, blue, and ma- 
genta represent gold, sulphur, hydrogen, carbon, nitrogen and iron, 
respectively. 



error the physical description will be qualitatively correct. We 
will use in this article the popular DFT+C/ approach, devel- 
oped by Anisimov and co-workers 16-18 , that represents a sim- 
ple mean-field way to correct for strong correlations in sys- 
tems with transition metals or rare earth compounds. 

The DFT+C/ method assumes that electrons can be split 
into two subsystems: delocalized electrons, that can be treated 
with traditional DFT, and localized electrons (3d or 4f), which 
must be handled using a generalized Hubbard model hamil- 
tonian with orbital-dependent local electron-electron interac- 
tions. The DFT+C/ functional is defined then as: 



£ DFT+C V(^), {n CT }] = E DFT [p a {~r^)] + 
+E mih [{n a }]- E DC [{n a }] 



(1) 



where i? DFT is the standard DFT functional; is the 

charge density for the a spin; _E Hub is the on-site coulomb 
correction; E BC is the double counting correction, which 
is necessary to avoid including again the average electron- 
electron interaction that is included in _E DFT ; and {n a } are 
the atomic orbital occupations corresponding to the orbitals 
that need to be corrected. 

The generalized Hubbard Hamiltonian is written as 



it. sr - - u — j N - . A 



(2) 

Following Ref. [l9|, we take the atomic limit of the above 
hamiltonian where the number of strongly correlated electrons 
N a = J2 m n mM is an integer and write: 



E DC = (integer N a \H- mt (integer N a ) 
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In contrast, for a noninteger occupation number, correspond- 
ing to an ion embedded in a larger system, we write: 



E = (noninteger N a \ H lnt | noninteger N t 
U ^ U-J 
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The above two equations can be merged and after some alge- 
bra the DFT+C/ functional can be written as 



2 ^ r 



• ^m,<r(i Tim, cr) (5) 



where U c s = U — J. 

It has to be noted that the correction term depends on the 
occupation number matrix. This occupation number matrix, a 
centrepiece of the DFT+C/ approach, is not well defined, be- 
cause the total density charge can not be broken down into 
simple atomic contributions. Since the appearance of the 
DFT+C/ approach, there have been many different definitions 
of the occupation number matrix 20 . We evaluate the occupa- 
tion number matrix by introducing projection operators in the 
following way2£: 
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(6) 



where ipj ^ are the KS eigenvectors for the j state with spin 
index a and qj are their occupations. The choice of the 

projection operators Pm ^ is crucial, because due to the non- 
orthogonal basis set different projection operators give differ- 
ent occupation number matrices. In our case, the choice cor- 
responds to so-called full representation. Here, the selected 
operator is 



P£ ] = \Xm)(Xm\ 



(7) 



where \xm) are the atomic orbitals of the strongly correlated 
electrons. Introducing this projector in (6), we get: 



full 



SD^S 



(8) 



where S is the overlap matrix and D CT is the density matrix of 
the system. 



III. COMPUTATIONAL METHOD 

We have performed the electronic structure calculations us- 
ing the DFT code SIESTA 21 , which uses norm conserving 
pseudopotentials and a basis set of pseudoatomic orbitals to 
span the valence states. For the exchange and correlation 
potential, we have used both the local density approxima- 
tion (LDA), as parameterized by Ceperley and Adler 22 , and 
the generalized gradient approximation(GGA), as parameter- 
ized by Perdew, Burke and Ernzernhof 23 . SIESTA param- 
eterizes the pseudopotentials according to the Troullier and 
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TABLE I: Energy gaps of metalloporphyrins in the gas phase, calcu- 
lated with LDA and GGA and with or without U. 





LDA 


LDA+U 


GGA 


GGA+U 


FeP 


0.10 


0.85 


0.48 


1.10 


CoP 


0.27 


1.80 


0.90 


1.80 


NiP 


1.55 


1.65 


1.50 


1.60 


CuP 


0.85 


1.40 


0.75 


1.40 


ZnP 


1.60 


1.60 


1.60 


1.60 



Martins prescription and factorizes them following Kleyn- 
man and Bylander 25 . We included non-linear core correction 
in the transition metal pseudopotentials to correctly take into 
account the overlap between the valence and the core states 26 . 
We also used small non-linear core corrections in all the other 
elements to get rid of the small peak in the pseudopotential 
close to the nucleus when using the GGA approximation. 

We placed explicitly the s and d orbitals of gold as valence 
orbitals and employed a single-^ basis (SZ) to describe them. 
We used a double-^ polarized basis (DZP) for all the other 
elements (H, C, O, N, S and TM). We computed the density, 
Hamiltonian, overlap and the matrix elements in a real space 
grid defined with an energy cutoff of 400 Ry. We used a single 
fc-point when performing the structural relaxations and trans- 
port calculations, which is enough to relax the coordinates and 
correctly compute the transmission around the Fermi level in 
the case of gold electrodes. We relaxed the coordinates un- 
til all forces were smaller than 0.01 eV/A. We varied the U- 
parameter and the radii of the projectors of the [/-projectors 
and compared the results with experiments or previous simu- 
lations of the isolated molecules. In this way we found that the 
optimal values for them were 4.5 eV and 5.5 A, respectively. 

Fig. (fl} shows the central part of the extended molecule in 
a gold-Fe-Porphirine-gold junction. The gold electrodes were 
grown in the (001) direction, which we took as the z axis. The 
sulphur atoms were contacted to the gold surfaces in the hol- 
low position at a distance of 1.8 A. We carried out a study of 
the most stable position and distance and found, both in GGA 
and LDA, that the hollow configuration was more stable than 
the top and bottom, in agreement with previous results ob- 
tained with other molecules. The most stable distances were 
1.6 A/ and 1.8 A/ for LDA and GGA, respectively. We finally 
chose a distance of 1.8 Afor all cases in order to make a sys- 
tematic study of geometrically identical systems. 

We performed the transport calculations using GOLLUM, 
a newly developed and efficient code 27 . The junctions were 
divided in three parts: left and right leads and extended 
molecule (EM). This EM block included the central part of 
the junction (molecule attached to the gold surfaces) and also 
some layers of the gold leads to make sure that the elec- 
tronic structure was converged to the bulk electronic structure. 
The same general parameters as in the bulk simulation (real- 
space grid, perpendicular A; -points, temperature, etc.) and also 
the same parameters for the gold electrodes (bulk coordinates 
along x and y and basis set) were used in the EM simulation. 



TABLE II: Magnetic moments in units of /xb of metalloporphyrins 
in the gas phase and between gold leads, calculated with LDA and 
GGA and with or without U. 





LDA 


LDA+U 


GGA 


GGA+U 


FeP 


2.00 


2.00 


2.00 


2.00 


CoP 


1.00 


1.00 


1.00 


1.00 


NiP 


0.00 


2.00 


0.00 


2.00 


CuP 


1.00 


1.00 


1.00 


1.00 


ZnP 


0.00 


0.00 


0.00 


0.00 


AuFePAu 


1.14 


1.32 


1.07 


1.34 


AuCoPAu 


0.62 


1.07 


0.83 


1.08 


AuNiPAu 


0.00 


0.00 


0.00 


0.00 


AuCuPAu 


0.84 


0.97 


0.80 


0.98 


AuZnPAu 


0.02 


0.00 


0.02 


0.00 



IV. ELECTRONIC STRUCTURE OF 
METALLO PORPHYRINS IN VACUUM 

Firstly we studied the effect of the U term on the electronic 
structure of the metallo-porphyrine MP molecules in vacuum. 
From these initial simulations it is already possible to fetch an 
idea of the main effects of correlations and their future impact 
on the transport properties. We chose for these analyzes a cu- 
bic unit cell of size 35 x 15 x 35 A to ensure that the molecules 
did not interact with adjacent images. 

In order to determine the influence of different ITs and cut- 
off radii and see how the results compare to previous exper- 
iments and calculations we studied first the case of the bulk 
oxide FeO 28 , where DFT is known to give qualitatively differ- 
ent results (metallic instead of insulating character) 29 . We per- 
formed calculations with U = 4 eV and U = 4.5 eV, which 
is the range of values most used in the literature for iron 29 , 
and we used projectors with different radii. We found that 
the results were very sensitive to these radii. Specifically, the 
gap for FeO only appeared for radii larger than 2.5 Bohr. The 
parameters that best fitted the experiments and previous cal- 
culations for FeO were U = 4.5 eV and r c — 5.5 Bohr, which 
gave a gap of 2.5 eV (2.4 eV in Ref. HH. 

To reproduce previous theoretical results published for the 
iron porphyrine 6,30 , we had to choose r cut = 1.5 Bohr, which 
is smaller than the values used for FeO. This is possibly jus- 
tified by the fact that the values of U and the projectors radii 
depend on the environment where the strongly correlated atom 
is located. In case of the other transition metals (Co, Ni, Cu 
and Zn), U is expected to increase towards the Zn but still be 
similar 17 . We therefore used the same parameters for all met- 
alloporphyrins, which also simplified the comparison between 
different cases and made the study more systematic. 

We summarize in the following subsections the main fea- 
tures of the electronic structure of all molecules. In some 
cases we show the density of states projected (PDOS) onto the 
d orbitals and/or the carbon and nitrogen atoms to determine 
the properties of the d states and their relation to the other 
molecular states. The gaps and magnetic moments of each 
molecule calculated with LDA and GGA and with or without 
U are summarized in Tables HI and ITR respectively. 
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FIG. 2: (Color online) Projected density of states (PDOS) on the 



(d) and d x 2_ y 2 (e) 



iron d states — d xy (a), d yz (b), d 3z 2_ r 2 (c), 
— for a iron metalloporphyrin in the gas phase, computed with GGA 
and spin polarization. Left (1) and right (2) columns correspond to 
calculations with U — eV and U — 4.5 eV. Positive and negative 
values represent spin up and spin down electrons. 



FIG. 3: (Color online) Projected density of states (PDOS) on the iron 
(a), nitrogen (b) and carbon (c) states for a iron metalloporphyrin in 
the gas phase, computed with GGA and spin polarization. Left (1) 
and right (2) columns correspond to calculations with U = eV and 
U = 4.5 eV. Positive and negative values represent spin up and spin 
down electrons. Notice the vertical scale is different in the middle 
panels (nitrogen). 



A. Iron porphyrine FeP 

The lowest energy electronic configuration for FeP in 
all calculated cases (LDA, LDA+U, GGA, and GGA+U) 
is [■■■]{d xy ) 2 (d z 2) 2 (dy Z ) 1 (d x 2_ y 2) 1 . This configuration is 
a consequence of the crystalline field generated by the 
molecule, which is located on the xz plane -the largest interac- 
tion between the nitrogens and the d states corresponds to the 
d xz , which moves up in energy and empties-. The total spin is 
therefore S z — 1, an intermediate-spin configuration, which 
is in agreement with experimental results 31 . The PDOS on 
the iron d states calculated with GGA and GGA+U is shown 
in Fig. (O for spin up and spin down electrons. As can be 
seen the closest orbitals to the Fermi level are the spin down 
d xy and d yz . The states that are filled and contribute to the 
magnetic moment of the molecule are the d yz and the d x 2_ y 2, 
whereas the d xz is completely empty. The d yz shows also 
a strong hybridization. When the U is included all gaps in- 
crease as a consequence of the movement of the filled orbitals 
to lower energies and the empty orbitals to higher energies. 
From here it is already possible to get an idea on the effect 
of the iron states on the transport properties of the molecule. 
States with large hybridization with other molecular states, 
such as the d yz are expected to give rise to relatively broad 
transmission resonances, whereas those states with small hy- 
bridization such as the (d z 2, and the d x 2_ y 2) are expected to 
produce either very thin resonances or Fano resonances, as 
explained below. 

The amount of hybridization can also be seen by plotting 
the PDOS on each type of molecular atom, as shown in Fig. 
(0 for GGA. Notice there are iron states that do not hybridize 
with the rest of the molecule, whereas other iron states do 



hybridize producing extended molecular orbitals. Focusing 
on particular states it is possible to see that the HOMO and 
LUMO have Fe and C contributions. In the HOMO, the con- 
tribution of Fe is bigger but in the LUMO both Fe and C con- 
tribute equally. Specifically, the Fe contribution to the HOMO 
(LUMO) comes from 3d xy and (3d yz ) orbitals; in addition, 
iron contributes to the HOMO-1 with the 3d z 2 and a small part 
of the 3d x 2_ y 2. In the case of carbon the state that contributes 
to both the HOMO and LUMO is the 2p y . Note that there 
are also sulphur states around the Fermi level, specially on the 
LUMO and below the HOMO but in order to simplify the de- 
scription we focus only on the carbon, nitrogen and metallic 
states. In this case the HOMO-LUMO gap is about 0.5 eV. 
With GGA+U the states that most contribute to HOMO are 
still Fe 3d xy , but the C 2p y states are the most important in 
the LUMO; the LUMO has also a small contribution from N 
2p y . Even though the U term acts only on the iron atom, the C 
2p y states are indirectly affected, i.e. the C 2p y states, located 
around -1 eV and 1 eV, lower their energy as a consequence of 
the changes on the iron states. In this case the HOMO-LUMO 
gap increases to 1.1 eV. 

With LDA and LDA+U the results are slightly different. 
The states that contribute most to the HOMO and LUMO are 
Fe M Z 2 and C 2p y . The HOMO-LUMO gap without U is 
about 0.1 eV. With LDA+U the Fe 3d z 2 states contributes 
much more to the HOMO and the C 2p to the LUMO, and 
the C 2p states are less affected by the U. The gap with U 
increases to 0.8 eV. 

By using the spatial projection of the density of states (lo- 
cal density of states, LDOS) it is also possible to understand 
where a particular molecular state is located on the molecule. 
In Fig. © we show the LDOS projected on the molecular 
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FIG. 5: (Color online) Projected density of states (PDOS) on the 
nickel d states —d xy (a), d yz (b), d z i (c), d^z (d) and d x 2_ y 2 (e) — 
for a nickel metalloporphyrin in the gas phase, computed with GGA 
and spin polarization. Left (1) and right (2) columns correspond to 
calculations with U — eV and U — 4.5 eV. Positive and negative 
values represent spin up and spin down electrons. 



plained in next section. 

FIG. 4: (Color online) Spatial distribution of the density of states 
projected on the spin-down main peak of the iron d yz (a) and d xz 
(b), computed with GGA. 



states associated to the d yz and d xz spin down states (with an 
isosurface value of 0.001 e/A 3 ), i.e. those where the weight of 
these d orbitals for spin down is largest, calculated with GGA. 
These states are located at E - E F — 0.29 eV (yz) and 2.29 
eV (xz) and move across the Fermi level when the atomic 
charge of the metal increases towards the Zn. Notice that in 
the case of the d yz state, due to the stronger hybridization, the 
d peak splits in two and therefore the choice is ambiguous. 
The spatial distribution of each peak is similar however. On 
each d state it is possible to see the typical shape associated to 
it, i.e. four lobes on the diagonals of the Y Z or XZ planes, 
plus some charge on other atoms of the molecule due to hy- 
bridization with other molecular orbitals. Notice that the d yz 
state does not interact too much with the nitrogens but rather 
with the carbons, specially with those located closest to the 
sulphur atoms. This produces an hybridization between this 
orbital and the carbon tt states, as can be clearly seen in Fig. 
© where the charge on the carbon atoms is mainly located 
on top of them. On the other hand, in the d xz case, the lobes 
are directed exactly towards the nitrogen atoms and therefore 
this state interacts strongly with them. This interaction is a- 
like, which is the type of interaction that nitrogen forms with 
the adjacent atoms, and therefore localizes the charge between 
atoms. The differences between the d yz and d xz orbitals and 
their coupling to different molecular states has a direct impact 
on the transport properties of some of these molecules, as ex- 



The lowest energy electronic configuration for CoP, in all 
cases, is [■■■](d xy ) 2 (d z 2) 2 (d yz ) 2 (d x 2_ y 2) 1 , which produces 
a ground state with S z = 1/2, in agreement with previous 
results' 1 . With GGA, the HOMO is made of Fe 3rf,,~ and 
3d xy states, C 2p y states and small fraction of N 2p y states. 
The LUMO is made only of Co 3d z 2 and a small fraction of 
3d x 2_ y 2. The gap is about 0.9 eV. With GGA+C/, in general, 
all states (Fe and C) lower their energy. Now the HOMO and 
LUMO are formed by C 2p y and N C 2p y states and have a 
small contribution of Co 3d xz states; the GGA +U is 1.8 eV. 
With LDA the results are similar, but the HOMO-LUMO gap 
is 0.3 eV, which increases to 1.7 eV for LDA+J7. 



C. NiP 

Unlike the FeP and CoP cases, the ground state elec- 
tronic configurations of NiP changes when the U is in- 
cluded. For GGA and LDA, the electronic configu- 
ration is [...](d xy ) 2 (d z 2) 2 (d yz ) 2 (d x 2_ y 2) 2 , i.e. S z = 

0, which makes the molecule diamagnetic. When the 
U is introduced, the electronic configuration becomes 

1. ^{d xy ) 2 (d z 2) 2 (d yz ) 2 {d x 2_ y 2) 1 (d xz ) 1 . The spin is S z = 1, 
and the molecule becomes magnetic. This is produced by the 
transfer of one electron from the 3d x 2_ y 2 to the 3d xz orbital. 
This can clearly be seen in Fig. ©, where without U both spin 
up and spin down 3d xz states are above the Fermi level (panel 
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FIG. 6: (Color online) Projected density of states (PDOS) on the 
nickel (a), nitrogen (b) and carbon (c) states for a nickel metallopor- 
phyrin in the gas phase, computed with GGA and spin polarization. 
Left (1) and right (2) columns correspond to calculations with U — 
eV and U = 4.5 eV. Positive and negative values represent spin up 
and spin down electrons. Notice the vertical scale is different in the 
middle panels (nitrogen). 



(dl)) whereas both 3d x 2_ y 2 are below the Fermi level (panel 
(el)). However, when U is introduced one of the 3d xz moves 
downwards and one of the 3d x 2_ y 2 moves upwards, each of 
them crossing the Fermi level. It can be therefore concluded 
that Hund's interaction is enhanced when U is introduced due 
to the larger repulsion between electrons. 

Fig.© shows the PDOS for Ni, C and N. With GGA the 
HOMO is formed by Ni, C and N states. Specifically, by Ni 
3d z 2 and 3d x 2_ y 2 and by C and N 2p y states. Unlike previous 
cases, the states below the Fermi level with spin up compen- 
sate the states with spin dow and the molecule becomes dia- 
magnetic. The LUMO, is formed by Ni 3d xz states and has a 
small contribution of N 2p x and 2p z states. The gap is about 
1.55 eV. However, when U is introduced, the HOMO is just 
made of C and N 2p y states. The LUMO is made of Ni 3d z 2 
and 3d x 2_ y 2 and by C and N 2p y states, just like the HOMO 
with GGA. In addition, as said before, the molecule becomes 
magnetic with S z = 1. The gap with GGA+U is similar to 
the case without U (1.65 eV). 

Similar words can be said when comparing LDA with 
LDA+C/. With LDA the molecule is again diamagnetic, with 
a gap of 1 .55 eV, like GGA. In this case however, the HOMO 
is only formed by Ni 3d states and in the LUMO the C 2p y 
states have stronger weight. With LDA+C/, the HOMO is 
formed by C and N 2p y states, and the LUMO by Ni 3d z 2 and 
3d x 2_ y 2 and by C and N 2p y states, just like with GGA+C/. 
The gap in this case does not change (1.55 eV). 



TABLE III: Change in Mulliken populations, Ae, when the metal- 
loporphyrins are coupled between gold leads, calculated with LDA 
and GGA and with or without U. 







LDA 


LDA+C/ 


GGA 


GGA+C/ 




AuFePAu 


0.097 


0.194 


0.134 


0.211 




AuCoPAu 


0.115 


0.260 


0.191 


0.286 


Spin up 


AuNiPAu 


0.255 


-0.729 


0.281 


-0.708 




AuCuPAu 


0.164 


0.236 


0.181 


0.277 




AuZnPAu 


0.273 


0.276 


0.301 


0.297 




AuFePAu 


0.545 


0.360 


0.529 


0.380 




AuCoPAu 


0.495 


0.246 


0.401 


0.263 


Spin down 


AuNiPAu 


0.254 


1.269 


0.276 


1.290 




AuCuPAu 


0.444 


0.315 


0.459 


0.309 




AuZnPAu 


0.273 


0.276 


0.297 


0.296 



D. CuP 

In CuP, the lowest energy electronic con- 
figuration is the same with and without U : 
[...'\{d xy ) 2 {d z 2) 2 {d yz )' 2 {d x 2_ y 2) 2 {d xz ) 1 - b . In this case, 
the 3d orbitals have closed shells, except the 3d xz , which 
loses a bit of charge. The total spin of the molecule is 
S z = 1/2, which mainly comes from the Cu 4s and 3d xz 
orbitals. Furthermore, the magnetism is not strictly localized 
on the Cu, but extends a little to the nitrogens. 

The PDOS on Cu, C and N with GGA, shows that the 
HOMO comes from C and N 2p y states. The Cu states are 
present in the HOMO-1, which is made of Cu 3d xz and has 
small contributions from N and C 2p x and 2p z states. The 
LUMO is made of Cu 3d. xz and N 2p x states. The energy gap 
is 0.75 eV. With GGA+C/, the HOMO remains the same, but 
the previous HOMO-1 moves to lower energies and the new 
HOMO-1 has C and N 2p y character. Although the LUMO 
is still formed by the same type of orbitals (i.e. Cu 3d xz and 
N 2p x ) now changes its spin polarization and becomes pop- 
ulated by spin down electrons, unlike GGA. The energy gap 
with U is about 1 .4 eV. 

With LDA the PDOS is similar to the GGA case, since the 
nearest states to Fermi level are the spin up Cu 3d xz , N and 
C 2p x and 2p z states. We therefore consider the HOMO to be 
located, as before, in the carbons and nitrogens. The HOMO- 
LUMO gap between spin-down states is 0.85 eV, which is now 
smaller than the HOMO-LUMO gap between spin-up states, 
1 .40 eV. The LUMO is the same as with GGA and is populated 
again with spin-down electrons. The results with LDA+C/ are 
the same as with GGA+C/, but shifted slightly in energy. 



E. ZnP 

The lowest energy electronic structure for ZnP is a close- 
shell electronic configuration (all the Zn 3d and 4s states are 
filled) due to the fact that the Zn states are all very deep in en- 
ergy as a consequence of the large nuclear attraction. There- 
fore the U correction does not affect to electronic configura- 
tion of this molecule. By analyzing the PDOS for Zn, C and 
N it can be shown that the HOMO and LUMO come from N 
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FIG. 7: (Color online) Projected density of states (PDOS) on the iron 
d states — d xy (a), d yz (b), d z 2 (c), d xz (d) and d x z_ y 2 (e) — for a 
iron metalloporphyrin between gold electrodes, computed with GGA 
and spin polarization. Left (1) and right (2) columns correspond to 
calculations with U = eV and U — 4.5 eV. Positive and negative 
values represent spin up and spin down electrons. 



FIG. 8: (Color online) Projected density of states (PDOS) on the 
nickel d states — d xy (a), d yz (b), d z 2 (c), d xz (d) and d x 2_ y 2 (e) 
— for a nickel metalloporphyrin between gold electrodes, computed 
with GGA and spin polarization. Left (1) and right (2) columns cor- 
respond to calculations with U = eV and U = 4.5 eV. Positive 
and negative values represent spin up and spin down electrons. 



and C 2p y states. The energy gaps in all cases are about 1.65 
eV. The only difference between LDA (LDA+C/) and GGA 
(GGA+C7) is the energy difference between levels, which is 
slightly different. 



V. ELECTRONIC STRUCTURE AND TRANSPORT 
PROPERTIES OF METALLOPORPHYRINS BETWEEN 
GOLD ELECTRODES 

A. Electronic structure 

When the metalloporphyrine molecules are coupled to the 
electrodes the most important effect on the PDOS is the broad- 
ening of molecular orbitals into resonances, as a consequence 
of the coupling and hybridization of the molecular states with 
the gold states. Although this effect is not very large on the 
d states it can still be seen by comparing Fig. (f2]i and Fig. 
(01. In case of iron, the largest difference is seen in the d xy 
PDOS (panel (b)), which spreads and decreases its height 
much more. This seems to indicate a larger derealization of 
this state as a consequence of its coupling to other molecular 
states that hybridize with the gold states. Something similar 
happens to the d xz and, to a lesser extent, to the other d or- 
bitals. This effect is maintained when the U is included. Also, 
as a consequence of hybridization and charge transfer the total 
spin of the molecule is reduced to almost a half of the value 
of the isolated molecule. This reduction comes mainly from 
the d xy state, which spreads and crosses the Fermi level. The 
charge transfer is summarized in Table UTT1 which shows the 
change in the Mulliken populations of the molecule when it 
is coupled between gold leads. As can be seen, all molecules 



gain charge, excluding the 'pathological' case of Ni. Since the 
gain is larger for spin up than for spin down the molecule ends 
up decreasing its magnetic moment, which is consistent with 
more delocalized states that reduce Hund's interaction. 

In metalloporphyrins other than FeP, the broadening of res- 
onances is not so clear, which can be explained by taking 
into account that the d states become more localized when 
the atomic number increases. In some cases the d states seem 
to even become more localized as can be seen by comparing 
panels (a) and (b) in Fig. © and Fig. dHJ. The case of Ni is 
also the most striking since the differences between the cou- 
pled and the isolated molecule are not only quantitative but 
qualitative. As can be seen in Fig. (JHJ, there is no splitting 
between up and down levels both with and without U, as op- 
posed to the case with U in the isolated molecule. This can 
be explained by arguing that electrons in the d xz and d x 2_ y 2 
can be now more delocalized, which decreases the effect of 
Hund's interaction. As a consequence, the total magnetic mo- 
ment with and without U in this case is 0.00 /is. as can also 
be seen in Table HJ 



B. Transport properties 

The zero-bias transport properties of these molecules are 
shown in Figs. ((5), ( fTOb and (fTTI) . calculated with LDA and 
spin polarization, GGA and spin polarization and GGA with- 
out spin polarization, respectively. As can be seen, increas- 
ing the atomic number from Fe to Zn produces qualitative 
changes in those transport properties at the beginning of the 
series. The main difference between the different metallopor- 
phyrins is due to two types of resonances that move to lower 
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FIG. 9: (Color online) Transmission coefficients for Fe (a), Co (b), 
Ni (c), Cu (d) and Zn (e) metalloporphyrins between gold electrodes, 
computed with LDA and spin polarization. Left (1) and right (2) 
columns correspond to calculations with U = eV and U — 4.5 
eV. Continuous and dashed lines represent spin up and spin down 
electrons. 



energies as the atomic number increases. One is a Fano-like 
resonance, a typical interference effector—, which is shown as 
a sharp increase followed by a steep decrease of the transmis- 
sion. The other seems to be a sharp Breit-Wigner resonance. 
These resonances appear in all cases, specially for Fe and Co, 
with GGA or LDA, with or without U and with or without 
spin polarization. The only differences are their position and 
width. With LDA and spin polarization, Fig. j9j, the Fano- 
like resonance appears clearly with Fe and, much sharper and 
on lower energies, with Co. This resonance moves to higher 
(lower) energies when U is added to Fe (Co). Above this res- 
onance there is a Breit-Wigner-like resonance that moves to 
lower energies from Fe to Zn and appears just at the Fermi 
level on CuP. It moves also to higher or lower energies when 
the U is included. With LDA and spin polarization, Fig. ( fTUb . 
the results are similar. Without spin polarization, Fig. (fTTT i 
we show for simplicity just the GGA cases. Again the behav- 
ior is similar. Some of these last results without U are also 
analogous to those obtained by Wang et al£L 

Obviously the two types of resonances must come from the 
d orbitals, due to their evolution with the metallic atom and 
also the fact that they are not present in porphyrins without 
metallic atoms. One hint about their origin can be obtained 
by looking at the projected density of states of the d levels 
and the surrounding atoms where these resonances happen. 
This is shown in Fig. (TTZt for the case of the Cobalt met- 
alloporphyrin, computed with LDA and without U, which is 
the configuration where different contributions from different 
atoms cam be most clearly seen. First, the spin down d states 
that are closest to the energy of the resonance are the d xz and 
d yz in the Breit-Wigner-like and Fano-like cases, respectively. 
For the first resonance the PDOS shows that there is a rather 
strong hybridization between the d xz level and the nitrogen 
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FIG. 10: (Color online) Transmission coefficients for Fe (a), Co (b), 
Ni (c), Cu (d) and Zn (e) metalloporphyrins between gold electrodes, 
computed with GGA and spin polarization. Left (1) and right (2) 
columns correspond to calculations with U — eV and U = 4.5 
eV. Continuous and dashed lines represent spin up and spin down 
electrons. 
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FIG. 11: (Color online) Transmission coefficients for Fe (a), Co (b), 
Ni (c), Cu (d) and Zn (e) metalloporphyrins between gold electrodes, 
computed with GGA and without spin polarization. Left (1) and right 
(2) columns correspond to calculations with U — eV and U = 4.5 
eV. 



atoms. For the Fano-like resonance, however, the d yz hy- 
bridizes much more with the carbon atoms. This is in agree- 
ment with the LDOS of Fig. which shows that the d xz 
hybridizes with the nitrogens and the d yz with the carbons. 

Taking into account the previous data it is possible to ex- 
plain the behavior of these systems as follows. The d orbitals 
are localized states that couple to certain molecular states. 
Such configuration is similar to that where a side state couples 
to a molecular backbone, which produces Fano resonances in 
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FIG. 12: (Color online) Transmission coefficients (a) and PDOS (b) 
around a Breit-Wigner-like resonance (1) and a Fano resonance (2), 
calculated on cobalt metalloporphyrin between gold electrodes and 
computed with LDA and spin polarization. Positive and negative 
values in the lower panels represent spin up and spin down electrons. 
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the transmission coefficients. The d orbitals generate there- 
fore Fano resonances associated to one or various molecular 
states. If, for example, the molecular state is the LUMO and 
the d orbital is in the HOMO-LUMO gap, the effect of the 
Fano resonance is clearly seen because it affects the transmis- 
sion in the gap, which has a large weight from the LUMO. The 
Fano resonance does not go to zero, however, because of the 
tails of the transmission of other molecular states. If, on the 
other hand, the d orbital couples to a state below the HOMO 
or above the LUMO, the effect of the Fano resonance turns 
out to be much smaller because the transmission of such state 
does not affect the transmission in the HOMO-LUMO gap. 
There is however one effect produced by the Fano resonance, 
which is related to its peak. This peak has a transmission of 1 
at the tip and therefore can be seen above the background of 
the transmission of other states. We can therefore confirm that 
the sharp Breit-Wigner-like resonances that are seen at differ- 
ent energies come in reality from Fano resonances. In view of 
these results it is possible to conclude then that the Fano fea- 
tures inside the HOMO-LUMO gap are produced by d orbitals 
that couple to the tt states, appear mainly on the LUMO. One 
d orbital responsible for such a feature is clearly the d xy , as 
seen in Fig. On the other hand, the sharp resonances are 
generated by d orbitals, such as the d xz , that couple to sigma 
states above the LUMO or below the HOMO. 

These features, specially the sharp resonance, appear in 
all cases and are a general characteristic of these molecules. 
We propose that the above-discussed metalloporphyrine 
molecules could be used as atomic or gas sensors. This is 
so due to the sharp changes in the zero-bias conductance 
that they generate when the Fano resonances cross the Fermi 
level, which should be altered whenever a given atom or 
molecules couples either to the metallic atom or to other parts 
of the molecular backbone. These resonances induce also dif- 
ferences between the spin-up and spin-down transmissions, 



FIG. 13: (Color online) Schematic representation of the model used 
to reproduce the ab-initio results. Case (a) corresponds to the d level 
coupled to the a state (labeled in the text as Hi) and case (b) to the 
d level coupled to one of the n levels (-H2). 



which are specially important in the case of iron, cobalt and 
copper. These could give rise to spin-filtering properties. The 
inclusion of the U term reduces however the spin-filtering be- 
havior, which only remains in the case of FeP. 



C. Simple model 



The previous behavior can be described with a simple 
model that takes into account the coupling of the molecular 
orbitals placed in the neighborhood of the Fermi energy to fea- 
tureless leads displaying a flat density of states, as well as the 
coupling of one of those molecular orbitals to the d orbital re- 
sponsible for the Fano resonanc o 38 i 39 . Specifically, the model 
includes four molecular levels: a a level below the HOMO, 
a level associated to the two linker sulphur atoms that repre- 
sents the HOMO and a n level which represents the LUMO. 
Finally, the d level has no direct coupling to the electrodes, but 
is instead coupled either to the a or to the it levels. The model 
and its two possible couplings are sketched in Fig. ( fT3l > (a) 
and (b), respectively. The Hamiltonian is therefore diagonal, 
except for the coupling between a given state and the d level. 
We also choose diagonal and identical V matrices, that couple 
the molecule to the leads: 
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To facilitate the comparison with the ab-initio results, we have 
chosen the following values for the on-site energies and the 
couplings: e a = -2, e s = -I, e d = 0, = 0.6, T a = 0.04, 
Ts = 0.06, T d = and r x = 0.06. where all energies 
are measured in eV. For model (a), we choose t a = —0.4 and 
t w = 0, while for model (b) we choose t a = and £ w = —0.2. 
The transmission and the retarded Green's function are then 
given by 



T(E) = Tr TG R \E)TG R (E) 



G R (E) 



EI - H aib - ^^ 



(10) 



We show in Fig. ( TBI the transmissions obtained for mod- 
els (a) and (b). The dashed line in Fig. (TBl i (a) shows the 
transmission coefficient found when the sulfur and it orbitals 
in model (a) are dropped and only the d and a orbitals are con- 
sidered. A clear-cut Fano resonance emerges at the d level on- 
site energy. However, this resonance is masked when the the 
sulfur and it orbitals are re-integrated back into the calculation 
leaving what looks at first sight a conventional Breit-Wigner 
resonance, the dashed line in Fig. (fT4l (b) is the transmis- 
sion obtained in model (b) when the a and sulfur orbitals are 
left aside, which features again a Fano resonance. However, 
the Fano resonance is now clearly visible in the full model, 
because the position of the a and sulfur resonance combined 
with the ordering of the divergencies in the Fano resonance 
can not mask completely the drop in transmission at higher en- 
ergies. The similarity with the ab-initio results for the cobalt 
porphyrine shown in Fig. (O is striking, despite the simplicity 
of the model. 




E-E F (e\) 



FIG. 14: (Color online) (a) Transmission coefficients calculated for 
model (a) described in the text and in Fig. d!3t (a). The black solid 
line corresponds to the full model (a), whereas the red dashed cor- 
responds to a simplified (a) model, where the S and n orbitals have 
been dropped from the calculation; (b) Transmission coefficients cal- 
culated for model (b) described in the text and in Fig. d!3t (b). The 
black solid line corresponds to the full model (b), whereas the red 
dashed corresponds to a simplified (b) model, where the S and a 
orbitals have been dropped from the calculation; 



electrodes. We have found two types of features in the trans- 
port properties that we show to be Fano resonances by the use 
of a simple model. We propose the use of metalloporphyrine 
molecules as possible nanoelectronic devices with sensing and 
spin-filtering functionalities. 



VI. CONCLUSIONS 

We have studied the electronic properties of isolated met- 
alloporphyrins and determined the influence of the exchange- 
correlation potential and strong correlations. We found that 
the HOMO-LUMO gap is greatly improved when a LDA+C/ 
or GGA+J7 description is used. However, we have found 
that the spin of the molecule does not change, excluding the 
case of NiP. We have also studied the electronic and transport 
properties of metalloporphyrins between gold electrodes. We 
found that the coupling to the electrodes changes only slightly 
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